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Abstract. We consider a Human Immunodeficiency Virus (HIV) model with a logistic 
growth term and continue the analysis of the previous article [6]. We now take the viral 
diffusion in a two-dimensional environment. The model consists of two ODEs for the 
concentrations of the target T cells, the infected cells, and a parabolic PDE for the virus 
particles. We study the stability of the uninfected and infected equilibria, the occurrence 
of Hopf bifurcation and the stability of the periodic solutions. 



Over the past thirty years, there has been much research in the mathematical modeling of 
Human Immunodeficiency Virus (HIV), the virus which causes AIDS (Acquired Immune 
Deficiency Syndrome). The research directions have been twofold: (i) the epidemiology of 
AIDS; (ii) the immunology of HIV as a pathogen. We are interested in the latter approach. 

The major target of HIV infection is a class of lymphocytes, or white blood cells, known 
as CD4 + T cells. When the CD4 + T-cell count, which is normally around 1000 mm 3 , 
reaches 200 mirT 3 or below in an HIV-infected patient, then that person is classified as 
having AIDS. 

Mathematical models have been proved valuable in understanding the in vivo dynamics 
of the virus. A gamut of models have been developed to describe the immune system, its 
interaction with HIV, and the decline in CD4 + T cells. They have contributed significantly 
to the understanding of HIV basic biology. 

Recently, the effect of spatial diffusion has been taken in account in HIV modeling. 
Funk et al. J7] introduced a discrete model: they adopted a two-dimensional square grid 
with 21x21 sites and assumed that the virus can move to the eight nearest neighboring 
sites. K. Wang et al. [18| generalized Funk's model. They assumed that the hepatocytes 
can not move under normal conditions and neglected their mobility, whereas virions can 
move freely and their motion follows a Fickian diffusion. In (TJ, two of the authors con- 
sidered a two-dimensional heterogenous environment: the basic reproductive ratio is gen- 
eralized as an eigenvalue of some Sturm-Liouville problem. Furthermore, in the case of an 
alternating structure of viral sources, the classical approach via ODE systems is justified 
via a homogenized limiting environment. 

In this article, we consider a HIV model which takes the viral diffusion into account 
in a homogeneous two-dimensional environment, and includes a quadratic logistic growth 
term as previously proposed in lfl6l[T7ll to consider the homeostatic process for the CD4 + 
T-cell count. The model reads: 
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1. Introduction 
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— = N/J,I - flyV + d v AV. (1.3) 

at 

The spatial domain is denoted by Q.( = (0, €) x (0, €), periodic boundary conditions are 
prescribed for V . Since the system ( ll.lt - (TT3l defines a dynamical system or semiflow, we 
will also use the abstract notation X(f) = (T(t), 7(f), V(t)). 

Our aim is to continue the analysis of the previous paper [6 1, where we studied the sys- 
tem jl.lt - dl.3t when dy = 0. We refer to [6| for an extended introduction to the biological 
issues. In brief, we recall that T and I denote the respective concentrations of uninfected 
and infected CD4 + T cells. The concentration of free virus particles, or virions, is V (for 
the sake of simplicity, we call V the virus). In dl.U . r is the average specific T-cell growth 
rate obtained in the absence of population limitation. The term 1 - T/T max shuts off T-cell 
growth as the population level T max is approached from below. Here fij is the natural death 
rate of CD4 + T cells, the term yVT models the rate at which free virus infects a CD4 + T 
cell. The infected cells die at a rate /// and produce free virus during their life-time at a 
rate N. In addition, fiy is the death rate of the virus. According to the literature (see e.g., 
TableQ]or [2j where \ij — 0.01, fij = 0.39), we assume the following biologically relevant 
hypothesis: 

Hj>Ht. (1.4) 
Note that the quantity r - the net T-cell proliferation rate, needs not to be positive (see 
lUSlP- 86]). 



Table 1 . Parameters and Variables 





Parameters & Variables 


Values 


Dependent 






variables 






T 


Uninfected CD4 + T-cell population 


mm -3 


I 


Infected CD4 + T-cell density 


mm -3 


V 


HIV population size 


mm -3 


Parameters & 






Constants 






r 


Proliferation rate of the CDA + T-cell population 


0.2 day' 1 


N 


Number of virus produced by infected cells 


1000 


a 


Production rate for uninfected CDA + T cells 


1.5 day~ l mm~ 3 


7 


Infection rate of uninfected CD4 + T cells 


0.001 day^mm 3 




Maximal population level of CD4 + T cells at 
which the CD4 + T-cell proliferation shuts off 


1500 mm- 3 




Death rate of uninfected CD4 + T-cell population 


0.1 day- 1 


/'/ 


Death rate of infected CD4 + T-cell population 


0.5 day- 1 


Hv 


Clearance rate of free virus 


10 day' 1 


Derived 






variable 






To 


CD4 + T-cell population for HIV negative persons 


mmT 3 



From a mathematical viewpoint: 

(i) N > and r > are parameters; 

(ii) the quantities (with associated dimension) a, y, fij, fir and fiy are fixed positive num- 
bers throughout the paper; 

(iii) r max is a large perturbation parameter, larger than any finite combination of a, y, ///, jxj 
and fiy of the same dimension (mm- 3 ). In particular, this hypothesis contains the condition 
T max > a/fi T of [161 P- 85]. 

The paper is organized as follows: in Section|2l we prove that System ( ll.lt -( fT31 > admits, 
for any value of the parameters r and N, the uninfected steady state X„ = (T u , 0, 0) and that, 
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in a region of the space of parameters, there exists also another steady-state solution, the 
so-called infected steady state X, = (T,, V,), where T,, /, and V, are positive. In the 
parameter space, we define the regions IX and 3 (this latter being the region where the 
infected steady-state exists), respectively for uninfected (the reproductive ratio is such that 
Ro < 1) and infected (Rq > 1). We recall that the basic reproductive ratio denotes the 
average number of infected T cells derived from one infected T cell (14)). We prove that 
the uninfected steady state is asymptotically stable in U, and unstable in 3. 

In El, we have exhibited an unbounded subdomain 7 in 3 in which the positive infected 
equilibrium becomes unstable whereas it is asymptotically stable in the rest of 3. In this 
unstable region, the levels of the various cell types and virus particles oscillate, rather 
than converging to steady values. This subdomain 7 may be biologically interpreted as a 
perturbation of the infection by a specific or unspecific immune response against HIV. In 
Section[3] we consider the linearization around X, = (T,, /,-, V,) of System ( ll.ll )-( fT31 l with 
Jacobian matrix A modal expansion of the resolvent equation enables us to construct a 
finite number of subdomains 7% (k = 0, . . . , Kq) in 3, that form a monotone non-increasing 
sequence (for the inclusion) with CPo = V. It turns out that the infected equilibrium X, = 
(Tf, If, Vi) is asymptotically stable for (N, r) e 3 \ 7 and unstable in the interior of 7. 
Therefore the stability issue is governed by the 0-th mode, hence similar to the case without 
viral diffusion (dy - 0). As a matter of fact, we are unable to confirm Funk et al. Q, who 
suggested that the presence of a spatial structure enhances population stability with respect 
to non-spatial models (see also 0]]). 

In Section H] we take the logistic parameter r as bifurcation parameter and prove the 
existence of Hopf bifurcations at the boundary 87. Since the system is only partially 
dissipative, the resolvent operator associated to the realization L, of ££i is not compact and 
therefore the proof demands more attention: it relies on the analyticity of the semigroup 
exp(fL,) (see e.g., IflOl . lfT31 ). Next, we perform a nonlinear analysis at the Hopf points 
via the Center Manifold theorem. It turns out that the bifurcating periodic solutions are 
independent of the space variables. 

Numerical illustrations are presented in Section[5] Finally, for the sake of completeness, 
we recall in an Appendix some basic facts about the eigenvalues of the two-dimensional 
Laplace operator with periodic boundary conditions and some Sturm-Liouville operators. 

Notation. For any i > we denote by L 2 the usual space of square-integrable functions 
/ : (0, if -> R. The square (0, if will be simply denoted by Q f . By (k = 1,2,...) 

we denote the closure in H k (the subset of L 2 of all the functions whose distributional 
derivatives up to k-th order are in I?) of the space of all k-th continuously differentiable 

functions / : R 2 — > R which are periodic with period I in each variable. The space H\ is 
endowed with its Euclidean norm. If X is any of the previous spaces, we write Xq to denote 
the space of complex-valued functions / such that Re / and Imf are in X. The norm in Xq 
is defined in the natural way: ||/H| = ||Re/H| + ||Im/H|. If v is a vector of C 3 , we denote 
by v\, V2 and V3 its components. Similarly, if f is a function defined in Q.( with values in 
R 3 (resp. C 3 ), we denote by f\, fj and fj its components. If the components of the vector 
v are complex numbers, we denote by v the vector whose components are the conjugates 
of the components of v. The Euclidean inner product in L 2 £ x L 2 e x L 2 e is denoted by (•, O2, 



i.e., 




for any f , g e L 2 e x L 2 e x L 2 e . Finally, we denote by Id the identity operator, and by (-) + the 
positive part of the number in brackets. 
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2. Equilibria 

In this section we are devoted to determine the non-negative equilibria of System ( II. U - 
( TOl l. i.e., the solutions (T, I, V)eL 2 xL 2 x Hi to the system 



max 



a-fi T T + rT[l-— \-yVT=0, (2.1) 

IT 



yVT-fi I I = 0, (2.2) 

NhjI-HvV + d v AV = 0. (2.3) 

To state the first main result of this section, let us introduce some functions and a few 
notation. 

By X„ and X, we denote, respectively, the function whose entries T u , I u and V„ are given 

by 

T u = T , i u = o, V„ = 0, 

where 



r-n T + ■sjir-H T ) 2 + ^ 



To = T (r) = T ma 

2r 

and the function whose entries are given by T,-, /, and V,-, where 

T = — 

'' yJV' 

j _ a fI T fIy | fi v r I Hv 



III y^]N y^N \ yNT, 



f*v 7 7\ jNT mm 
We further introduce two sets which will play a fundamental role in all our analysis, 
namely the uninfected and infected regions U and 3 in the parameter space, which are 
defined by 

U = \(N, r) : N > 0, r > 0, R (N, r) < 1}, 3 = \(N, r) : N > 0, r > 0, R (N, r) > 1}, 
where 

Ro(N, r) = 

is the reproduction ratio. 

The interface Ro{N, r) — 1 between the two regions It and 3 is the graph of the mapping 

Hv L T ~r+ J (fi T ~ r) 1 + ) 

AUW = — i \ = -^r, (2.4) 

2ay 7io('") 

which is decreasing by virtue of the condition r max > afi T l . Its image is the interval 

{^—, Inverting the roles of and r, it is useful to define the inverse mapping: 

rait(N) = -, — r , N > N crit (+oo) - 



As it has been stressed in the introduction, throughout the paper we assume that 

Mi>Mt- (2.5) 
To prove the following theorem we assume also that 

T lmx >T^ x (a,y,fij,iJ T ,iJ V ), (2.6) 
where T max is fixed and large, and depends only on the quantities in brackets. 
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Theorem 2.1. The following properties are satisfied: 

(i) in U the uninfected steady state X„ = (T U ,I U , V u ) is the only non-negative equilib- 
rium; 

(ii) in 3 there exist two non-negative equilibria, respectively X„ and the infected steady 
state X ; = (Ti,Ij,Vi); 

(iii) it holds < T u < T max and, in 3, < T, < T mm ; 

(iv) X„ and X, are all the possible equilibria to System ( ll.U - (11.3l l, with all the compo- 
nents being non-negative. 

Proof (i) Suppose that X = (T, I, V) is a solution to System (l2~TTi-(l231l. Then, from (l2~TT i 
we deduce that 



■Mr-yV+ ^j(r-fj. T -yV) 2 + Ti 



4ar 

max 



T = -V : ,,v. (2.7) 

2r 

Replacing the expression of / given by ( 12.2b in (12.31 > and, then, using (12.71 ). we obtain the 
following self-contained nonlinear equation for V: 

d v AV - p v V + <D(V) = 0, (2.8) 

where 

yNV 



O(V) 



2r 



4ar 

■p T -yV + J(r-p T - yV) z + 



T 

1 max J 



As it is easily seen any solution to ( 12.81 ) in H 2 leads to a solution to System (11.1 b - (fT73T >. 
Moreover, from any non-negative solution to Equation ( 12.81 ) we can obtain an equilibrium 
to System ( ll.ll )-( fT31 ) will all the components non-negative in D,{. Hence, we can limit 
ourselves to looking for non-negative solutions V e //? to Equation ( 12.8b . 

Clearly, ( 12.81 ) admits the trivial function V = as a solution. This solution leads to the 
equilibrium X„. 

(ii) Let us look for other positive constant solutions to Equation ( 12.81 ). We are thus lead 
to look for solutions to the equation <E>(V) - fiyV = which are non-negative. 

A straightforward computation reveals that V, is the unique solution to such an equation. 
Moreover, for any fixed r > 0, V, is positive if and only if > N CI i t (r) (see ( 12.41) ) i.e., if 
and only if (r, N) e 3. In this case, replacing V = V, into ( 12.71 ) and ( 12.21 ), we immediately 
conclude that the function X, is an equilibrium of System dl.ll )-( fT3] ). 

(iii) Showing that T u < r max is just an exercise. On the hand, the inequality T, < r max 
in 3 follows from the definition of T, observing that, if (r, AO e 3, then 



yT u yT max 

(iv) To prove that X„ and X, are the only equilibria of System ( ll.U -( 11.3l ) with all the 
components being non-negative, we adapt to our situation a method due to H.B. Keller 
fl3l . We argue by contradiction. We suppose that X = (T, I, V) is a solution to System 
( 12.U - d23b with V non-negative in £2f and not identically vanishing, and such that V + V,. 
Let us set W := V — Vj. Since both V and V, are solutions to ( 12.8b . clearly W e and 
solves the equation 

d v AW - fi V W + A(V, Vi)W = 0, (2.9) 

where 



yN 



Mx,y) = 

2r 



r-fi T -yx+ J(r-p T - yx) 2 + 



4<3T 



y 2 NT, 



max 



T 

* max J 

y 2 NT m . dx f 1 y(tx + (l - t)y) - r + p T 



2r 



-y 



J[r-p T -y(tx + (l-t)y)f + ^ 

y 1 max 



for any jc, y > 0. 
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Let A max (V, V/) and A lmx ( V, 0) denote the maximum eigenvalues in L 2 of the operators 
d v A + A(V, Vj)Id and d v A + A(V, 0)ld, respectively. By Corollarv lA.21 we know that 



A n 



s (^0) = - inf 



and 



We now observe that 
A(V, V,)- A(V,0) : 



inf 



'd v j ni \Vif/\ 2 dxdy - j nt A(V, 0)if/ 2 dxdy 
' dv Li Wtyfdxdy ~ f n[ MY, Vi)ip 2 dxdy 

L t ^ 2dxd y 



y 2 NT n 

2r 



, y 2 NT n 

2 2r 



-V, 



-V, 



-1 + 



-1 + 



f 

Jo 



y 2 NT n 

2r 



■V; 



1 



y(fV + (l -t)Vi)-r + fi T 
^[r-ii T -v{tV + {\-t)V i )] 2 + ^ 

\y(tV + (l-t)Vi)-r+fi T \ 
y j[r-ti T -y(tV + (l-t)V i )] 2 + ^ 

yiWWo. + vo + r + fiT 



zdt 



dt 



c. 



A /[r + ^ + r(IIVIU + y i )] 2 + ^^ 

since the function x i— > x(x 2 + 4ar/7' max )~ 1 '' 2 is increasing in [0, +oo). Here, we have taken 
advantage of the Sobolev embedding theorem to infer that V e H 2 is continuous in flf . 
Note that the constant C is positive. From this remark we can easily infer that 

^max(V,V/) </Wx(V,0)-C. 

Since W satisfies ( 12.91 ) and it does not identically vanish in Qf, //y < A max (V, V,) < 
/W(V,0). 

To get to a contradiction, we now rewrite the equation satisfied by V in the following 
way: 

d v AV + A(V, 0)V - A max (V, 0)V = Qi v - WV, 0))V := Z. 

Fredholm alternative implies that Z should be orthogonal to the function iff which spans the 
eigenspace associated with the eigenvalue A m \ n ( V, 0) of the operator dyA + A( V, Q)Id. But 
this can not be the case. Indeed, by Corollary I A. 21 the function if/ does not change sign in 
Of. Moreover, since V is non-negative in D.( and it does not identically vanish in Q.( and, in 
addition, fiy < A mix ( V, 0), Z is non-positive and it does not identically vanish in Q^. Hence, 
if/ is not orthogonal to Z. □ 



3. Stability of the equilibria 

In this section we are going to study the stability of the equilibria X„ and X,. We begin 
by studying the stability of the uninfected equilibrium X„. 

Theorem 3.1. The following properties are satisfied: 

(i) in the domain 11 the uninfected equilibrium X„ is asymptotically stable; 

(ii) in the domain J the uninfected equilibrium X„ is unstable. 

Proof. To avoid cumbersome notation, throughout the proof we do not stress explicitly the 
dependence of the functions and operators on r and N. 

We prove the statement showing that the linearized stability principle (see e.g., IfTTl 
Chpt. 5, Cor. 5.1.6]) applies to our situation. For this purpose, we begin by observing 
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3 


Ro > 1 








u 


R <1 





Figure 1. Profile of the curve r i-» N cr it( r ) (i- e - - 1) which defines 
the two domains U and 3. With the values of Table Q] N cr it decreases 
from (i T fiy/ay = 666.67 to nvlyT m!LX = 6.67. 



that, for any (N, r) the linearization around X„ of Problem (ll.ll l-( fT73l is associated with the 
linear operator £„ defined by 



((>■ 



A, = 



l*T - Y^t) Id 








-Hi Id 
wNId 



-jTJd 

yTJd 

dyls. - [l^/Id, 



Its realization L u in (L^) 3 with domain D(L U ) = I: x x Hj c generates an analytic 
strongly continuous semigroup. Indeed, L u is a bounded perturbation of the diagonal op- 
erator 

([r-UT -PA Id 

\ 1 max / 

-nild 

d v A - livid, 

which is clearly sectorial since all its entries are. Hence, we can 

■2 
tt,C 



A 



defined in I;, x /. • x H? c 

apply [ 14 Prop. 2.4. l(i)] and conclude that L u is sectorial. Since is dense in L£, the 
associated analytic semigroup is strongly continuous. 

To complete the proof, we need to study the spectrum of the operator L u . We fix A e C 
and consider the resolvent equation AX - L„X = F, where X = (T, I, V) e Lg x I? c x H? c 

and F = (F\,F<i,Fi) is a given function in l. : x /.\ x L^. Writing the previous equation 
componentwise gives 



[r-UT - p^)T -yT u V 
-mI + yT u V = AI-F 2 , 



AT - Fi 



(3.1) 



/i/JV/ + c/vAV - yUyV = AV- F 3 . 



If /I ^ - / u / we can use the second equation to write / in terms of V. Substituting it in the 
last equation we get the following self-contained equation for V: 



d v AV -\A + nv- 



A+H! 



V = -F 3 - 



A+fij 



F 2 . 



(3.2) 



x 
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We recall that the spectrum of the realization A of the Laplace operator in Lg, with 
as a domain consists of eigenvalues only and it is given by 

4 * 2 -a , a, 



<r(A) = ^-—(kt+fy:k u k 2 €: 

(see AppendixlAli. Hence, if c(A) :- A+py- 71 %+J" does not belong to cr(A), then Equation 
(13.2b admits a unique solution V e H% c - A straightforward computation shows that c(A) is 
real if and only if A is real. Moreover, the function A i-> c(A) is strictly increasing in [0, +00) 
and c(0) = fx v - yNT u is positive if (N, r) e U. Hence, if (N, r) e It, then c(A) ( cr(A) for 
any A with non-negative real part, and Equation (13.2b is uniquely solvable. 

We can now uniquely determine / € L^, from the second equation in ( 13.11 ). Finally, from 
the first equation in ( 13.11 ). observing that 



2rT u / , 4or 

r-fir- ~ — = - J(r - ht) 2 + - — < 0, 

* max V * max 

we can uniquely determine T e L?,. We have so proved that any A e C with non-negative 
real part is in the resolvent set of the operator L„, if (N, r) € U. In view of the linearized 
stability principle this implies that the trivial uninfected solution to System ( U.U -l TOI ) is 
asymptotically stable. 

Let us now suppose that (N, r) e 3 and prove that L„ admits an eigenvalue with positive 
real part. As above, we are led to consider the function A h-> c(A). Now, c(0) < 0. 
Hence, there exists A* > such that c(A*) = € <x(A). This implies that Problem ( 13.11 ). 
with F\ — F2 = F3 = and A = A*, admits a non trivial solution, i.e., A* € cr(L„). 
Again, the linearized stability principle implies that (T„, 0, 0) is unstable. This completes 
the proof. □ 

The issue of the stability of the infected solution is obviously much more complicated. 
For notational convenience we sort the eigenvalues of the realization A of the Laplacian 
in L 2 C , with H^ c as a domain (see Appendix lAll. into a non-increasing sequence {— Afr}. 
Similarly, we denote by the eigenfunction associated with the eigenvalue A^. This allows 
ut to expand any function /eLL into the Fourier series 



/=2> 



k=0 

where f% denotes the £-th Fourier coefficient (with respect to the system (e^)) of /. As it 
is observed in the proof of Theorem lA.il Aq is simple and all the other eigenvalues are 
semisimple, and their multiplicity can be computed explicitly. 
We are going to prove the following results. 

Theorem 3.2. Under the hypothesis ( 12.61 ), it holds: 

(i) the infected equilibrium X,- = (Tj, V,-) is asymptotically stable for (N, r) 6 5 \ T, 
where 7 is defined in ( 13.171 ) with k — 0; 

(ii) the infected equilibrium X; = (T;, V,) is unstable in the interior of 7. 



3.1. The resolvent equation. As in the proof of Theorem l3.11 we do not stress explicitly 
the dependence on r and of the operators and the sets that we consider in what follows. 
For (N, r) e J, the linearization around X, = (T,, V,) of System dl.lMl.3l l is associated 
with the linear operator 



fflld -fid 

l„ ^ + r(l-^)l« -mid fid 



t*yr 
'NT m 

Npild d v A - fiyld, 



(3.3) 
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Proposition 3.3. For any (N, r) e J the realization L, of the operator in (L^,) 3 with 
domain D(Li) = L 2 * X /.\ X ffj c generates an analytic strongly continuous semigroup. 
Moreover, the spectrum of Li is given by 

pyr ayN 



yNT mdx p v 
where, for any k e N, o~ k is the spectrum of the matrix 



km 



fyr _ ayN r, _/iy ^ 



N 



Afy// -fiM* - /iy. 



M t = 



Proof. The same arguments as in the proof of Theorem 13.11 show that L,- generates an 
analytic strongly continuous semigroup in {L^ff . 

Let us determine its spectrum. For this purpose we use the discrete Fourier transform. 
If a function v = (vi, v%, V3) in I:, x L 2 x H? c solves the resolvent equation Ay - L,v = f, 
for some A e C and f = (/i , /j, yj) in (L^) 3 , then its Fourier coefficients v* = (vijt, v 2,k, v Xk ) 
(k = 0, 1, . . .) solve the infinitely many equations (Aid - M k )y k = fit (k = 0, 1, . . .), where 
fit = (fi,kr fz,k< fi,k) an d /),* denotes the fc-th Fourier coefficient of the function f (j = 
1,2, 3). Clearly, any eigenvalue of M k (k = 0, 1, . . .) is an eigenvalue of L,. Therefore, 

Cr(Li) D UkeNO~k- 

On the other hand, if A £ cr k for any k = 0, 1, . . ., then all the coefficients (v\^, Vajfc, V3,i) 
are uniquely determined through the formulae 

vi,* = ^ 1 {[(/* + + <M* + y"v) - Wvl/i,* - Wvh,k ~ ^rO* + M/)/3,ife| , (3.5) 



1 f A + flfy/lj; + jUy 9 2 2 

V2,k = — L(-/V + a 7 N T max - yp T pvNT mdx + yp v rNT m . dx )f uk 

V k (A)( yp v NT mdX 



+ {Ayp v NT mdx + p\r + ay 2 N 2 T mdx )f ltk \ 



IA (2p v r - yrNT mdx + yp T NT mdx + AyNT m!ix )f Xk \, (3.6) 



yN 2 T, 

1 f A*/ 



v Xk = rn ttt I ~ — [(ffy"N r max - yp T pvNT mdx + yp v rNT m . dx - p v r)f\,k 

T> k (A)[yp v T mdx 

+ (Ayp v NT mdx + p 2 v r + ay 2 N 2 T mdx )f 2 , k ] 
+ -(Ayp v NT m . dx +p 2 v r + ay 2 N 2 T mdx )f 3 ,X (3.7) 



ypvNTn 
where 

T> k {A) = A 3 + d u A 2 + d 2 , k A + d 3Jc , k = 0,l,..., (3.8) 

and 

pyr ayN 

di, k = d hk (N, r) = d v A k + p r + p v + — - — + , (3.9) 

yNT m . dx py 

ayN Pvf 

d.2,k = d 2 , k (N, r) = pid v A k + ayN + (yU/ + d v A k ) + — — — (ja + p v + d v A k ), (3.10) 

P-v yNT m . dX 

d\ k = d Xk {N,r) = pip v {r - p T ) + ° WlN {p v + d v A k ) + -^^—{d v A k - p v )- (3.11) 

P-v yNT nvdX 

Note that, if A differs from both — 22^1 an( j then 

T> 1 r\ a \ ^ A ^ ^ _Eyr\ aypiN pip v r \ 
V k {A) ~d v \A + A \p, + + — — — + + — — — \ A k , 

{ \ p v yNT m . dX p v yNT mdX ) 
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as k — > +00. Hence, for any A £ {—■^f- ^7, u UteN °"£> ^ holds that 

Vl ' k ~ A 2 + AU + 2* + + 22El! + jw /U ' 

\r' fly jNT m ) fly 



V2,k 



V3,k 



(-r^ + ^^-^)fu + (^t + ^ + A)f 2 , 
1 

d v A k {A 2 + A^ l + ^ + ^) + ^ + ^-} 

2 ,, ,, ,. v /™„„, \t2 



x 7 — + - ViVtN / u + + — — + )f 2 ,k 

+ (A + AI ,)( A + 2^ + _ff!2_k 4 
V //y yNT mRX r 

as A: — » +00. Thus, the sequences fvi^}, {v2,jt} and {/Uvj^} are square-summable. This shows 
that the series whose Fourier coefficients are Vi^, v 2 ,k and v^jt, respectively, converge in I 2, 
(the first two ones) and in TfrL (the latter one). The inclusion cr c j — U 

v ' (t,C v ' l yNlmn fly ' ^ J J 

U&eN °"ifc follows. On the other hand, the previous computations show that, if ^ = -/// or 

/I = then the series having v;^, V2,* and as Fourier coefficients do not, in 

general, converge in L 2 (the first two ones) and in H 2 (the latter one). Hence, these values 
of A belong to the essential spectrum of L, . Thus, ( 13.3b is proved. □ 

3.2. Study of cr k . Clearly, at fixed k = 0, 1, . . ., each set cr k consists of at most three 
eigenvalues vj lk , j = 1, 2, 3, either all real, or one real and two complex conjugates, which 
verify the equation A 3 + d\ k A 2 + d 2 ^A + d^^ = 0, with dj ik (j = 1,2,3) being given by 

The Routh-Hurwitz criterion enables us to determine whether the elements of cr k have 
negative real parts. The latter holds if and only if d\ tk , ^3,* and the leading Hurwitz de- 
terminant D 2k = d\ tk d 2yk — d^ t k are positive. The case k — corresponds to the system of 
ODEs considered in j6j- 

Obviously, di )k > 0. As far as d^ is concerned, we remark that d3 )k (N, r) > d-3 t o(N, r), 
which is positive in J and vanishes at N a -i t (r). 

For (N, r) e 3, we compute the Hurwitz determinant D 2 ^(N, r) and we get 

° 2 ' k(N > r) = 2 lL T 2 ( A ^ 2 + B ^ N > + C kW) • ( 3 - 12 > 



max 



where 

A k = fiyfjxj + n v + d v A k ), 

Bk(N) = yiu 2 v NT mRX [fU V dlA 2 k + 2ayd v NA k + 2^ v d v A k + 2ynyL V d v A k 

- WifivNT m . ix + 2ayfi,N + 2ayfi v N + fxjfiy + 3/^i^ 2 v + /4]> 

C k (N) = N 2 y 2 Ti^(fi If i 2 v d 2 v A 2 k + a 7f i v d 2 v NA 2 k + a 2 y 2 d v N 2 A k + 2a 7f i 2 v d v NA k 

+ fiifiydvA k + i?niydyA k + 2ayniHvdvNA k + ayfiyN + ay/ui/jyN 
+ ay/j. 2 /i v N + a 2 y 2 /d V N 2 + a 2 y 2 fijN 2 + fijfiTfiy). 
Both A k and C k (N) are positive, whereas B k (N) vanishes at 

N _ Hv(d 2 v A 2 k + 2(hj + \xy)d v A k + 3/ji/jv +fi 2 + //y) 
y(2ad v A k + 2a^i + 2a/x v - fJ-iMvTmax) 
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Figure 2. Profile of the curve No yk . 



In ( I3.13I ). as a function of A k , the denominator vanishes at 

_ HiHvT mdx _ Hi + jiy 
lady dy 

which is positive and generically does not meet any of the A k s, k > 1. There are two cases 
(see Fig. 13: 

(i) < A k < Ao, hence N 0Jl > 0. Then, B k (N) > if < N < N 0Jc ; 

(ii) A k > Ao, hence No, k < and in this case B k (N) is positive for any N > 0. 

The sign of the polynomial A k r 2 +B k (N)r+C k (N) is obviously related to the discriminant 

A k (N) = (B k (N)) 2 - 4A k C k (N) = y 2 ^ v N 2 T 2 mx (a k N 2 + b k N + c k ), 
which in turn is of the sign of a k N 2 + b k N + c k . The coefficients a k , b k and c k read: 
a k =y 2 HiT m!ix {LiiiivT m a X - 4ad v A k - 4a{p,i + jiiy)}, 

b k = - 2ynin v \[d 2 v A 2 k + 2(/i/ + Hv)d v A k + + 3^/jUy + MyJ^max - 4ad v A k - 4a(jX] + ft v )}, 
c k =n v {{d\Al - /u 2 ) 2 + 4nvd\A\ + 6/i v (Mi + V-v)d 2 Y A 2 k 

+ 4fX V \Mv + 3 Wv + Mi( 2 Mi ~ Hr)]d v A k 

+ 2/^u y (3jU/ - 2n T ) + 6yU /y Uy + /j,/j 2 v (11/j, - 4fi T ) + /j 4 v 

Let us examine the signs of these coefficients, 
(i) As a function of A k , the coefficient a k vanishes at 

- A 



A 2 



4ady dy 4ady 

Clearly, A2 is positive thanks to (12.6b and, as Ao, generically does not meet any of 
the A k \ for k > 0. Then, a k is positive if < A k < A2 and negative otherwise. 

(ii) b k < due the hypothesis (12.61 1. 

(iii) q > under the biologically relevant hypothesis fii > yij (see ( 12.51 )). 
Next we compute: 

6 k :=b\ - 4a k c k 

= 16yV/yMy"/ + + d v A k ) 

x [a(dlA 2 k - yL//) 2 T' m ax + 4afi v d 3 v T m:ix Al + (fij/j. v T max + 6a// v + 4a^ I n v )d 2 ,T m . dx A 2 k 
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+ jUy[/Z;(jU/ + pv)T 2 mx + 4a{p] + 2p,p v + p 2 , - p///r)7 , ma x + 4a 2 p I ]d v A k 

+ (i}(i T Hv T ZwL + a h4 + fyiMv + h l v + fyjfivfai ~ Mr) + M/Mv(5yU/ - 4// r )]T max 

+ 4a 2 p 2 jp v + 4a 2 p I p 2 ,'\. 

Again, thanks to the hypothesis pj > pr, S k is always positive. Therefore, the roots of 
a/cN 2 + b k N + c k = 0, namely the ones of A k (N) - 0, are: 



Ni* = 



-b k - 

2a k 



N 2 ,k = 



-b k + 

2dk 



(3.14) 



Remark 3.4. Since Bk(No ik ) = and A k (No,k), C k {No,k) are both positive, a k (No tk )) 2 + 
bkNot + Ck < 0. This provides us with a useful information regarding the position of No,k 
with respect to Ni^ and N 2ik , according to ( 13.14-b . i.e., Ni tk < Nq^ < N2, k if a k > and 
Noji < N2, k or No tk > Ni ik if a k < 0. In particular, No, k can meet neither N\jc nor N 2 ^ k - 

Now we are in a position to begin the discussion, depending on the position of A k . Let 
us distinguish three cases. 



x10 




Figure 3. Comparison of the three curves Nq^, Ni, k and A^jfc, as functions of A k . 



Case I: < A k < \ 2 . In this situation a k > 0, b k < and c k > and k ranges in a finite 
set of indexes. It is an extension of the case k — 0, see (5). It follows from Remark [3~4l 
that B k vanishes between Ni,k and N 2 ^ which are both positive. In particular, for T m . dx large 
enough (as we are assuming), Ni£ > pv /(y?max)- Indeed, 

,2 



a k 



+ b k 



+ c k 



2,2 ,2x2 , .,.2 ,3,3 , A ,._. .2 ,2,2 , -r,,3 ,2 ,2 , 3 . 



=H V {d^A% - pjY + 4p l v d v A{ + 4p I p%d M v A l k + f>p?yd L v A\ + Sp^dvAk 
+ 4p I p 2 v {p I - Hr)d v h + 4p 4 v d v A k + 2p]p 2 v + 2p]p\ + 2p]p 2 v {p I - p T ) 
+ 4pip% + 4p I p 3 v (p I - p T ) + p\ + o(1), 



as T„ 



¥oo. Again, since pj > pj, 



+ b k 



+ Ck > 0, 



if T max is large enough. It follows that either p v /(yT m!lx ) < JV U or pvlijTm^) > N 2l k- But 
as it is immediately seen, N 2 , k > Hv/(yTmax)- Indeed, 

d 2 w A\ + 2{pi + p v )d v A k +p} + 3pipv + f*l 



7(J-iT B 
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2y 2 j^ 1 fi v T m . dx 

3yUy 



as T max — > +00. Hence, fiv/(yT mdx ) < Nij as it has been claimed. 

We consider four subcases depending on the position of N with respect to Ni, k and N 2>k . 

(i) Assume N hk < N < N 2 , k - Then, A k (N) < 0. Since A k > 0, it follows that D 2 ,k(N, r) > 
for all r > r cl ; t (A0 (see ( 13.12t and recall that we are taking (N, r) from 3). 

(ii) If N < N lk , then B k (N) > 0. Hence, A k r 2 + B k (N)r + C k (N) > for any r > since 
A k , B k , C k > 0. It thus follows that D 2 ,k(N, r) > for all r > r clit (N). 

(iii) If AT > N 2 , k , then B k (N) < 0. Since A k (N) > 0, the equation A k r 2 +B k (N)r+C k (N) = 
admits the two real and positive roots: 



-B k (N) - ^/A k JN) , m ~B k (N) + ^fA k W) n 1V . 

n,*(A0 = — , r u (N) = — . (3.15) 

LA k LA k 

Observe that 

2C k 



-B k + a/AZ 



^—riniHTl^v + ayn^N + ay^^N + ayfinyN + a 2 y 2 fi v N 2 

+ a 2 y 2 /uiN 2 + a 2 y 2 d v N 2 A k + 2ay/j 2 / d v NA k + ^i^? v d v A k 
+ niii 2 v d\A 2 k + aypL V d\NA\ + fi 2 fild v A k + 2aynin v d v NAk) 

>/i T > = r ait (N), (3.16) 

as r max — > +00. Consequently, the Hurwitz determinant D2, k {N,r) is positive for 
r C nt(N) < r < riji(N) and r > r 2tk (N), it vanishes at r — r\^{N) and r = r 2tk (N), and is 
negative for r\ jc(N) < r < r 2i k(N). 
(iv) Assume N e {N\jc,N 2 j^. In such a case, A k (Nj, k ) = and the polynomial A k r 2 + 
Bk{Nj, k )r + C k (Nj tk ) has the double root r k (N^k) = -B k {N ^ k )/2A k . However, this 
solution makes sense only if B k (Nj tk ) < 0. Hence, only the case N = N 2 ,k is relevant, 
and we have D 2 , k (N 2 , k ,r) > for r > r cdt (N 2 , k ) except at r(N 2 , k ) = -B k (N 2 , k )/2A h 
where it vanishes. 

We are now in a position to define the subdomain T k of 3 by 

y k = {{N, r):N> N 2 , k , r uk (N) < r < r 2 , k (N)} , (3.17) 

see (EH and fllfl . and at N = N 2 , k , r uk (N 2 , k ) = r u (N 2 , k ) = -B k (N 2 , k )/2A k . In the 
domain 3, d\jt and dj^ are positive, and D 2>k (N, r) is positive except in 3\. More precisely, 
the Hurwitz determinant D 2>k (N, r) is negative in the interior of y k and it vanishes on the 
boundary of y k . 

Case II: A 2 < A k < A . Now, a k < 0, b k < 0,q > 0. Hence, _/V 2 ,a < O.A^* > 0. Since 
A^o^ > 0, it holds that A^t > A^* according to Remark [3~4l Moreover, as it is immediately 
seen, Nq^ > (y^max) Vv- There are two possibilities: 

(i) (N, r) e 3 satisfies A^ < N , k . Then, B k (N) > 0. Therefore, A k r 2 + B k (N)r + C k (N) > 
for any r > r CI n(N) since the coefficients are all positive. It thus follows that 
D 2 ,k(N,r)>0. 

(ii) (N, r) e 3 satisfies N > N 0Jc . Then, A^ > N uk and a k N 2 + b k N + c k < 0. There- 
fore A k (N) < and A k r 2 + B k (N)r + C k (N) has the sign of A k which is positive, so 
D 2 AN,r)>0. 
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Case III: A k > A . Here, N^t < and, therefore, B k (N) > for all N > 0. The conclusion 
is the same as in Case II (i). 

We summarize our results in the following proposition. 

Proposition 3.5. Denote by K2 the largest integer such that A^ < A2. Then, 

(i) for k = 0, . . . , K2, the Hurwitz determinant DzjciN, r) is, respectively, negative in the 
interior of the subdomain y k of J, positive in J \ T k , and it vanishes on the boundary 
ofy k ; 

(ii) for k = K2 + 1, K% + 2, . . , the Hurwitz determinant D2, k (N, r) is always positive in J. 

Remark 3.6. To give an idea, with the numerical values of Tableland t — 1, A2 = 1239.5 
and lies between /I97 = I16n 2 and Agg — 1287T 2 , therefore K2 = 97. 

To conclude this subsection we prove the following proposition which gives a much 
clearer picture of how the sets 7 k are ordered in the space of the parameters. 

Proposition 3.7. Let K2 be as in the statement of Proposition U3\ Then, the following set 
inclusions hold: 

V Kl £ 0V 2 -i £ ■ • • c y t c . . . c c ? . 

Proof. To begin with we claim that N2,k < N2,k+\ ( see ( 13-141 )) for any k = 0, . . . , K% — 1. To 
prove the claim we observe that 5k = 6(Ak), where 

6(x) =16y 2 /iipv(Mi + P*v + dyx) 

x \adyT max x A + 4afi V d 3 v T m:lx x 3 + (fi 2 jp v T max + 6a/4 - + 4ap I ^ v )d 2 ,T m . dx x 2 

+ Hvltfip-i + P-v)^^ + 4a(f/j + 2p,fi v +p 2 v - pip T )T mzx + 4a 2 n I ]d v x 

+ yU^UrjUy^max + a \-t4 + fylfXy + fJ-y + 4fl 2 p.y(jJ, - p T ) + fXif£y(5(J.i - 4/l r )]r max 

+ 4a 2 /xjpv + 4a 2 pifi v }. 
We compute the derivative of the function 5 and get 

5'(x) = 16yV/yUv{5a<4r max x 4 + 4ff(^/ + 5fj. v )d 3 v T max x 3 

+ 3(ju 2 p v T m . dx + iafijpv - 2a/i 2 + \Qafi 2 v )d v 

+ a{9p. I p 2 v + 3fi 2 p v - l^mp-v 5//y)r max + 2o? '(j, I (j,y]d v x 

+ pjfiviprfiv + yUy + j"/ + 2/i / /i l /)7 ,2 ax 

+ a(17fi 2 fil - 8fiifi T fi v + 8yU///v + 16^1 + 5/4 + f4 ~ 8tfp.TP-v)T m ax 

+ Sa 2 p 2 p.v + &a 2 iiip v }dv- 

Under hypothesis ( 12. 6t this function is positive and, consequently, A: 1 — > is non-decreasing. 

Similarly, ak = a(Ak) and bk = b(Ak), the functions a and b being strictly decreasing. 
Hence, the sequences [ak] and {bk} are non-increasing. Since a k > and bk < we now 
easily get the claim. 

To complete the proof of the inclusion 3\ + 1 c 7 k for any k = 0, . . . , Ki— 1 , we show that, 
for any N > N2M\ we nave r \,k(N) < r iyk+ i(N) < rz >k +\(N) < ^(A). These properties 
follow immediately from the definitions of r\ k{N) and r2, k {N) observing that < Aj < 
Aj+i, Bj(N) < Bj +l (N) < (since N > No,h for any h = 0, . . .,k + 1; recall that we are in 
the Case I(iii) where N\ >k and N2, k are both positive, and take Remark [3~4l into account) and 
< Cj < Cj+i for any j = 0, . . . , K 2 - 1. 

Finally, since Aq < A t ,N2,\ > A^.o- Consequently, Ti is properly contained in 3V □ 
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Figure 4. The sets 3>£. 



3.3. Proof of Theorem l3.2l The proof follows from Propositions ^ .5! l3Jl Routh-Hurwitz 
criterion and the linearized stability principle. 

(i) As Proposition I3.5I shows, for k = 0, . . . , K2 the leading Hurwitz determinants 
D2jc(N,r) (see ( 13. 121 )) are positive in 3 \ 7 k . On the other hand, if k > K 2 + 1, then 
D 2 ,k(N, r)>0 for any (N, r) e 3. 

By Proposition 13.71 it holds that 7% c 7q, k = 1,2, Hence, we conclude that 

Di,k{N,r) > for any k € N, if (N,r) e 3 \ 7q. Since the other two Hurwitz determi- 
nants are positive in the whole of 3, it follows from the Ruth-Hurwitz criterion that, if 
(N, r) e 3 \ 7, then all the element of UfeN °"a nave negative real part. Hence, Re <r(L,) < 
(see d3.4t ). It remains to invoke the linearized stability principle as in the proof of Theorem 
HU 

(ii) The instability of X„ can be deduced from H which deals with System dl.lMl.3t 
in the case when dy = and shows that, in this situation, the infected equilibrium X, is 
unstable. □ 



4. HOPF BIFURCATION AND INSTABILITY 

For fixed N > Owe take the logistic parameter r > r cl ; t (A0 as a bifurcation parameter. 

We recall that at fixed (Af, r) e 3, System dl.lt -l TOb has two equilibria: the uninfected 
trivial solution X U (N, r) and the infected, positive solution X,(iV, r). At X,(A^, r), the Ja- 
cobian matrix is £, = &i,N,r, see d3.3t . As we already observed in Proposition 13.31 the 
realization L;,2v,r of the operator in (i|) 3 with domain D(Lj^ r ) - //; x I: x 
generates an analytic strongly continuous semigroup that we denote by e ,Li - N - r . 

In this section we are interested in proving that Hopf bifurcation occurs on the boundary 
of the set 7 (i.e., at the points (N, r u0 (N)) and (N, r 2 fl(N)) with N > N 2 ,o, where n$(N), 
r2,o(N) and A^.o are given by (13. 14-b and d3.15t ) and in analyzing the stability of the bifur- 
cated periodic solutions. 

Note, that for T max large, 

n] + 3mn v +v 2 v + 2 ylm^THvijdi + n v ) 1 

^2,0 = 7f — + o(T max ). 

J^I 1 max 
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Hence, A^o is positive if r max is large enough, let us say, if T > T^ ax (pT,^i,l^v, a, y) > 
T^. dx (p T ,pi,p v ,a,y) (see (12.61 )). We assume hereafter that 

T^>T x m!LX (a,y,pj,pT,Hv,N). (4.1) 

Here, differently from the previous sections, to avoid confusion we stress explicitly the 
dependence of the operators, numbers and sets that we consider on r. We do not stress 
the dependence on N since in the following discussions only the parameter r varies, N is 
(arbitrarily) fixed. In particular, we simply write r\ and r-i instead of r\{N) and r2(N). 



Theorem 4.1. Let N > A^.o be fixed. Under the hypothesis ( 14. U , Hopf bifurcation occurs 

at the critical points r — rj, j — 1,2. More precisely, 
(i) for any (3 e (0, 1), there exist cq > and smooth functions rj,pj : (sq, bq) — * R and 
X kj : (-e , eo) -» C 1+/J (R, L 2 xL 2 xL 2 )n C^(R, L 2 xL 2 x H 2 ) (j = 1,2) such that 
pj(0) — 1, rj(0) — rj, Xy(0) = X,-(ry), Xt/(c) is not constant in time ifc^O and its 
period is 2npj(x)la>j, where 



aypiN pvrj 

u= x \ayN+ H — — — (jii+nv), j = 1,2. 

^ p v yNT mRX 

(ii) There exists <5 > such that, ifX € C 1+/5 (R, L 2 x L 2 x Lh n C^R, LjxLjx Hp is a 
periodic solution to System (1 1 . 1 b - d 1 - 3b (vv/zere r is replaced by 7) with period Inp/a), 
such that 

||X - X / (r j )|lc/<(R,LjxLjxH, 2 ) + IIX - X.i(rj)\\ cl+HKL 2 xL 2 xL p + \r- rj\ + |1 - p\ < 6 , 

for j — or j — 1, then there exist s € (—eo, sq) and to € R such that T — tj and 
X = Xy(e). 

Proof. We limit ourselves to considering the case when r — r\, the case r — r2 being 
completely similar. 

For r in some neighborhood of ri, we set u = X - X,(r), s = r — r\ and write System 
( ll.U - dl.3t at the infected equilibrium as 

^=nn,s), (4.2) 
dt 

where 

<x < ^ (Mv(s + n) , uyN\ p v n+s 2 
ii(u,*J = — I 1 u\ M3 Mj — yu\uy, 



?2(U, S) = 



-fj. T + (s + n) l- 



^(u, s) -P1NU2 + dvA-UT, - pyU3 

Clearly, by the Sobolev embedding theorem, 2F is a smooth function defined in L 2 xL 2 xH 2 . 
Note that the derivative iF u (0, 0) is the operator L, n in Proposition l3.3l More precisely, 

L <>. = [^-Mr + ri (l- 5^-)] Id -pad %Id 

Npild d v A - p v Id) 

By Proposition 13.31 operator L;. n is the generator of a strongly continuous analytic 
semigroup in (Li) 3 . 

Let us prove that <x(L, n ) consists of eigenvalues with negative real part and a pair of 
purely imaginary and conjugate eigenvalues /li(ri) and feiri), which are simple eigenval- 
ues and satisfy the transversality condition. Once checked, these properties will yield the 
assertion in view of [ 14, Thm. 9.3.3] (which deals with fully nonlinear problems but, of 
course, it applies also to the semilinear case). 
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Being rather long, we split the proof into four steps. 

Step 1. Here, we prove that cr(L, ri ) consists of eigenvalues with negative real part and 
a pair of purely imaginary and conjugate eigenvalues. For this purpose, we observe that, 
since CP^ is properly contained in 7q for any k — 1, . . . , K 2 (see Proposition 13. 71 , the pair 
(N, ri(N)) belongs to 3 \ 7 k for any k — 1, . . . , K 2 . Therefore, from the results in Subsection 
I3.2l and Proposition l3.3l it follows that cr^ ri is contained in the halfplane {A e C : Re A < 0). 

As far as <tq is concerned, Orlando formula (see e.g., (8] Chpt. XV]) shows that the 
Hurwitz determinant D 2 ^(r) (see ( 13.121 )) factorizes as follows: 

O 2 ,o(r) = -(Mr) + A 2 (r))(A 2 (r) + /l 3 (r))(/li(r) + A 3 (r)), 

where Ai(r), A 2 (r) and A^r) are the roots of the polynomial 

Vo, r (A) = J- 3 + d U) (r)A 2 + d 2fi (r)A + d 3 .oW 



(see ( 13.81 )) (i.e. the elements of cr ). The point (r\,N) lies on the boundary of TV Hence, 
Difi(f\) vanishes, i.e., 

(Ai(n) + A 2 (r 1 ))(A 2 (r ] ) + -l 3 (n))(Ai(n) + Un)) = 0. (4.3) 

Since the coefficients of Do, ri are real and positive, at least one of the three roots A\(r{), 
Mji), ^3(^1) (let us say A 3 (r\)) is real and negative and the other two roots are either both 
negative or they are complex and conjugate. From ( 14. 3t it follows that Ai(ri) and A 2 {r\) 
are purely imaginary and conjugate. 

Step 2. Let us prove that there exists a gap between cr(L,' ri ) \ {A\(r), A 2 (r)} and the 
imaginary axis. We have to consider the set cr> = crk, ri (k = 1,2,...) which consists of the 
roots of the third-order polynomial T>k, ri (see ( 13.81) ). Indeed, as we have already remarked, 
Ai(ri) is negative. 

Write A = p. - M. If A is a root of the polynomial T>k, r , , then p. = A + M is a root of the 
polynomial p k , ri (A) = ^ 3 + di,k{r\)^ 2 + d z ,k(r\)A + d X k(r\), where 

di,k(r\) = di, k (r{) - 3M, 

d 2 , k (n) = d 2 ,k(n) - 2Mdi, k (r\) + 3M 2 , 

d3, k (n) = dxk(n) - Md 2 ,k(n) + M 2 d uk (r{) - M\ 

As it is easily seen 

di,k(n) = d v A k + o(A k ), 



d\k(.r\) = 



aypiN pipvn I ayN p v n \ 2 
+ ^—^ \pi + + — \M + M 



P-v jNT mKX \ p v yNT, 



d v A k + o(A k ), 



di, k (ri)d 2 ,k(n) ~ dxkin) = (^Li + — + ~^p~ ~ 2M) A\ + o(A 2 k \ 

\ p V yNlmux I 

as k — > +00. Hence, if we take M satisfying the inequalities 

Pv yNT m . dX \ p v yNT mm j 
Pv yNT miiX 

then, for A: sufficiently large (say k> K3 > K 2 ), d\.k(r\), d^j^ri) and d\^(r\)dxk(r\)-di t k(r\) 
are all positive. Hence, Routh-Hurwitz criterion applies and shows that the roots of p k , ri 
have negative real part of any k > Kj. As a byproduct, {Jk>K 3 °~k,n c {A e C : Re A < —M). 

Since Uis^a:, °~k,n consists of finitely many eigenvalues with negative real part, up to 
replacing M with a smaller constant if needed, we can assume that (Ji</t °"*,ri c {A e C : 
ReA< -M}. 

Step 3. We now prove that the eigenvalues A\(r\) and A 2 {r\) are simple. 
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First, we prove that the resolvent operator R(A,L^ ri ) has a simple pole at Aj(j\) (j = 
1,2). We limit ourselves to proving this property for the eigenvalue A\(r{), since for the 
other one the proof is completely similar. 

From the proof of Proposition l3.3l we know that, for any A e p(L, >ri ) and any f e (Li) 3 , 



f +co +00 +00 \ 

R(A, L ;>r , )f = vi, k (A)e k , ^ v u (A)e k , ^ v 3 ,k(A)e, 

\k=0 k=0 k=0 



where Vj,k(A) (j = 1, 2, 3, k e N) are defined by Q3l)-(f377l) 
Observe that 



aydyN fi v nd v \ . ay/j,d v N p.i^i v r\d v 
d v A + fijdy H 1 A h h 



A k + d(A, r\), 



fi v yNT m . dX ) n v yNT„ 

where d(A, r\) is independent of k, and it is smooth in A. 

As it is immediately seen the coefficient in front of A k does not vanish at A — A\(r\). 
Hence, there exist a neighborhood U of A\(r\), toeN and a positive constant^ such that 
|f fcfif/QI ^ X^k for any k > ko and any A in U. From ( 13.5l )-( 13.7t we thus deduce that 

IvjytU)! ^ c (l/ul + Lfr*l + IA*D, 

for any k> ko, j - 1,2,3 and /let/. 

Since D^C/liCn)) ^ for any k ^ 0, the previous estimate can be extended to any 
A: > 1. Hence, 

/ +00 +00 +00 \ 



R(A,L Ul )f = (vi, (A),v2,oW), v 3 ,o(^)) + 



Vyt=l i=l fc=l / 

2\3< 



where the second term in the previous splitting defines a function with values in £((L£) ) 
which is bounded in t/. 

The singularity of R(-,Li n ) at /I = Ai(ri) is due to the first term of the splitting. The 
results in Step 1 show that A h-> Do,r,W has a simple zero at A — A\(r{). It thus follows 
at once that the function A \-* ((A - A\(r\))\>i 7 o(A), (A - A\(ri))v 2 fi(A), (A - M(fi))v3,o(A)) is 
bounded around A = A\{r\). 

Summing up, we have proved that the function A h-> (A — A\(r\))R(A, L^) is bounded 
around A = Ai(r\). Consequently, R(-, L, ri ) has a simple pole at A — A\(ri), so that, by lfT4l 
Prop. A.2.2] Ai(ri) is a semisimple eigenvalue of L, ,-, . 

To conclude that it is, actually, a simple eigenvalue, we have to show that the eigenspace 
associated with A 1 (r\ ) is one dimensional. This property follows from recalling that Dk, ri (A\ (r\ )) 
if k > 1. Hence, any eigenfunction associated with Ai(r\) is a constant. 

Step 4 . We now check the transversality condition. Observing that 

di,o( r i) = -0*1 (ri) + A 2 (r l ) + i 3 Oi)) = A 3 (n), 
diflin) = A x (r l )A 2 {r l ) + A l (r l )A 2 (r i ) + A l (r 2 )A 2 {r i ) = A(r x )A(r 2 ), 
from (13. 9t and (13.1 Ob we conclude that 



AM) = (-lylyN + + -^(^v), J = 1.2. 

and 

/I3O1) = -Hi-pv ~ —rpp • 

By IIT21 Chpt. 20] the function r \-* /li (r) is smooth in a neighborhood of n . Hence, 
differentiating the formula (A\(r)) 3 + difl(r)(A\(r)) 2 + d 2 ${r)A\(r) + d 3 fl(r) = 0, evaluating 
it at r = r\ and then taking the real part, we get 

d? n (n) - d\ () (ri)d 2 fl(ri) - d\, o(n)d' Jr{) 
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The sign of (j^Re/li) (n) is the sign of d\ (r l )-d[ li (r x )d 1 < A (ri)-d\ Xl (ri)d' 1Q (rx). A straight- 
forward computation shows that 

^3,o( r i) ~ d' lfi (n)d 2 ,o(ry) - di i0 (ri)d' 20 (n) 
= 2 I 2 [rWrA^T^ ~ 7 NT max (3/i/My + 2ay^N + 2ayp v N + fi 3 v + tffiv) 

y ^ * max 

Since 

n 1 —(miMtM 3 v + ayp 3 v N + aypip^N + ayrffivN + a 2 y 2 ^ v N 2 + a l y 1 p. I N 2 '), 



v 



as r max — > +oo (see ( 13.161 1). ^Re/li) (rO is positive if T max is sufficiently large, as we are 
assuming. Hence, the transversality condition is satisfied. This completes the proof. □ 



Proposition 4.2. The bifurcated periodic solutions provided by Theorem \4.1\ are indepen- 
dent of the spatial variables, i.e., they are the same bifurcated periodic solutions of the 
following system of ODE' s: 

— =a-p T T + rT 1-— )-yVT, (4.4) 

at \ r max ) 

j^yVT-^f (4.5) 
dV 

— = Nful - nyV. (4.6) 
at 

Proof. In J6] it has been proved that System (I4.41 >-( f4~6b exhibits a Hopf bifurcation at 
r - r j U - 1» 2). A branch of periodic solutions bifurcates from X;(ry) ( j = 1,2). Clearly, 
such solutions are space independent. Moreover, a statement analogous to Theorem l4.1t ii) 
holds for the Hopf bifurcation associated with Problem ( 14.4l )-( Pk6i >. see ifTOl Thm. II, p. 
16]. Therefore, up to replacing eo with a smaller value, if needed, we can infer that, for 
any s € (-eo, eo), Xj(s) coincides, up to a translation in the time variable, with one of the 
bifurcated periodic solutions in [6 Thm. 4.5]. This shows that any function Xj(e) is space 
independent. □ 

We can now prove the following theorem: 

Theorem 4.3. Suppose that r max > T ( ^) dx , where depends on a, y, pij, fir and fiy (see 
the proof). Then, the following properties are satisfied. 

(i) If N < N* (where N* is the first positive zero of the function 3~C in ( 14.131 l), then the 
periodic solution X^(e) is orbitally asymptotically stable with asymptotic phase. 

(ii) For any N > 0, the periodic solution X^(e) is orbitally asymptotically stable with 
asymptotic phase. 

Proof. The arguments in Henry's book IfTTII (see also in a more general situation) show 
that the stability of the bifurcated periodic solutions can be read on a Center Manifold. 
This allows to reduce our problem, which is set in a infinite dimensional Banach space, to 
a problem in a finite dimensional space. 

To obtain this finite dimensional problem, we first need to determine the spectral pro- 
jection associated to the eigenvalues -tOji and a>ji (j = 1,2). As a general fact, such a 
projection is the sum of the spectral projections Pj i+ , associated to the eigenvalue iojj, and 
Pj-, associated to the eigenvalue -ia)j. Since io)j and -ia)j are simple eigenvalues (see 
Theorem 14. U . there exists a unique projection on the eigenspace relative to iutj which 
commutes with L, . Similarly, there exists a unique projection of the eigenspace relative to 
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-ia>j which commutes with L,. Using these facts it is easy to check that 



for any v e (L 2 ^) 3 , where 



and 



<Pj 



'jNT m , : 



fly + 1(0 



N(sj - iojj) fijN 



*J 



111 - ILOj 

' /UN 



ayN 



ayN 



~Ht+ 1- 



yNT tr 



Kj = (<Pj,lffjh, 



for j = 1,2. In what follows we set <pj = (<Pj t \,(pi2, <fij^) and if/j - (iffj,\,ifrj,2, fitf)- 

As it is well known, Pj allows to split (L^) 3 into the direct sum of the two subspaces 
Pj((L 2 c ) 3 ) and (I-Pj)((L 2 c ) 3 ) where Pj((L 2 c ) 3 ) = {z<p + wlp : Z, w e Q and (/ - Pj)((L 2 c ) 3 ) = 



{u e (Lf c ) : (v, 1^)2 = (v, 1^)2 = 0). In particular, Pj maps (L ) into itself and allows us to 
split the space (L 2 ) 3 into the direct sum of the two subspaces Pj((L 2 ) 3 ) = {z<p + Ijp : z e C) 
and (/ - Pj)((L 2 ) 3 ) = {u e (L 2 ) 3 : (u, ^2 = 0). 
Let us rewrite Problem ( 14.21 > in the form 



du 

— =Li, r/ u + S 7 -(u,5), 



(4.7) 



where 



b/,i(u, s) = — — — mi - 



rj + s _ 



yNT ma , 



S 7 -, 2 (u, s) =s 1- 



T 



u 1 — yu\uy, 



S 3 ,/u, s) =0. 

Splitting Problem gjj along Pj((L 2 ) 3 ) and (/ - Pj)((L 2 ) 3 ), we see that any solution u e 
C'([0, a) x (L 2 ) 3 ) n C([0, a)xL 2 xL 2 x Hh to Problem (|4~7T >. defined in some time domain 
[0, a), can be identified with the pair of functions (z, w), with z(f) e C and w(f) e L 2 xL 2 xH 2 
for any t e [0, a), which solves the system 



dz _ 

— = iW/Z + 9 ;(z, Z, W, s), 

off 

t/w _ 

— = L,w + 3C.-(z, z, w, s), 
dt 



(4.8) 
(4.9) 



where 

5j(z, Z, W, s) =A£j 



' n -7-, (TT^i^i) 

mi(i - <A>,i) 7 



t\ 2 



+ 2kj 



y(l -^ ( i)Re(^,i)- 



(0 + ^l^ilVj.i 



zz 



+ Kj 



f 

J a. 



2(r; + s)<fjii//j\ 

y(l - ^j,i)((pj,iW3 + wi) — ■ — -wi 



yNT m!iX 



f + sl- 



max 

a 



dxdy 
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/ 

JClr 



7(1 - ^,l)(^,lW3 + Wi) W\ 



max 



dxdy 



"I 



11 

7(1 - lffj,l)W3 Wi - 



yNT ms} 



+ s 



w\dxdy; 



9Cj(z, z, w, s) -Sj(z, z, w, s) - 2Re (§j(z, z, w, s)^>), 



for = 1,2. Modulo the identification of Pj((L 2 ) 3 ) with the set {(z,z) : z e Q, the Center 
Manifold for System (I4.8b - d479l > is the graph of a smooth function Yy of the variable (z, z, s), 
defined in a neighborhood of zero with values in (/ - P ; )((L 2 ) 3 ). 

The equation to be analyzed, to understand the stability of the bifurcated solutions Xj(s), 
is therefore the following one: 



dz _ _ _ 

— = icojz + S j(z, z, Tj(z, z, s), s) =: g/z, z, s). 



(4.10) 



This ODE can be studied with classical methods (see e.g., IflOl Chpts. 1 & 2]). One needs 
to expand the nonlinearity gj around as 



gJ (z,z,0)= J] 



2<h+k<3 



U) 

°hk - , /i i4\ 
7777 Z/A + 0(\Z\ ). 

h\k\ 



The coefficients gi,k(s) are fundamental to determine the stability of the periodic solutions 
to ( 14.101 ). In fact, such solutions are stable if and only if Re c\(rj) < 0, where (see e.g., IflOl 
P- 90]) 



= 2^1^-21.^-^; 



,01,2 



1. 



Si) a 



+ - 



(./) 

2 °2i- 



To expand g around the origin, one first needs to expand the function T( , •, 0) around 
(0, 0). Since this function is smooth, we can expand it as 

Tj(z, z, 0) = a\ j, z + a 2 U) z + a 3 °V + a 4 (j) zz + a 5 °V + 0(|z| 3 ). 

Replacing u(f) = z(t)<p + z{i)lp + T(z(f),z(0) into System (I4.8b -d4-.9l>. expanding 

%j(z,z, w, 0) = 3C y -iz 2 + 2Xj#z + 0C~[f + 0{\z\\w\) + 0(\w\ 2 ), 

and observing that 

T w(o = — kz{t\W))z'{t) + —i( Z (t),zJtj)zvj, 

dt dz oz 

an asymptotic analysis reveals that 

Tj(z,z, 0) = z 2 (2icjj - /.,,■) ';JC/.| - 2zzL^pCj, 2 + z\-2iu>j - /.,,. ) '!K~7 + 0(\z\ 3 ). 
where 



It 



max 



(2^ 2 Re (W,i^,i) - l) - 2y^-,iRe (jc^O 



• mi 



Re(Kj(p j)2 i/fj,i) + m,i (! _ 2 ^ 2 Re(^; 2 )) 



2t l 



+ m,i 



Re(^i)-2rr^,iRe(*r / ) 
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' (rjpj\_ + ^ ( \ ^ {2rs u ^ ; | ^-, _ jj _ 2 ^2 Re ( ^ . l)Re 



2 ^2 fej^ilL + rRe . j)] Re (A ;y -,,^,,) + yRe (<p A1 ) (l - 2{ 2 Re ( W , 2 )) 

\ * max / 



|2 



2^ | + yRe Re (^,i) - 2y^Re (^ ;1 )Re (*,-) 

* max 



Note that, since JCj(z,z,w) e (/ - Pj)((L 2 ) 3 ) for any z e C and w e (7 - f j )((L 2 ) 3 ), 
3Cy ( i, 3C y ;2 and UCyi belong to (/ - Pj)((L^) 3 ) and the operators ±ia> - L^ r . are invertible on 

(/-P,)('(Z2) 3 ). 

Now a long but straightforward computation shows that 



= t + 2 yKj<PjMia - tj.iY 



2 



max 



„(/>_ 2 - _ 



g (j) = _ -jr^ -jrj,- + 2 7K ^ {iff . 2 _ ^ . x)( \ 

1 max 

= " 2rjK / j ' l{2 [(ni(^j ' LirjT'OC^h - %l(^Ki2)l] 

-* max 

+ y*/^ - ^~7)f 2 [^-T((2/w j - /.,, ) + «2iajj - L U r 1 T x 'Kj,{)i 

where (•)* denotes the k-th component of the vector in brackets. 

Since an explicit computation of these coefficients for any value of T max is uneasy, and 
we are interested in large (enough) values of r max , as in J6] Sec. 4.3] we determine the sign 
of Reci(r,) via an asymptotic analysis as T m . dx — > +oo. We get 



Rec 1 (r 1 ) = ^5C(iV) + (l) > Re c x (r 2 ) = - 50 ^^ v) + (r m 2 x ), (4.11) 

max 



3 

!)(»)""" v ~" — 
where D(A0 and 5t(JV) are respectively given by 

D(A0 =2a(p.in v + fj 2 v + ayN)(fij fj 2 v + 2fjrfj\ + 6ayfiifj v N + fly + 6ayfi 2 v N + a 2 y 2 N 2 ) 
x (a 2 y 2 N 2 + 3ayfj.jfj.vN + 3ayfi 2 v N + fp\fj 2 v + 2fJifj\ + fj 4 v ) 
x (n]n v + 2fj If j 2 v + fiifj v + fj\ +fJ 2 v + ayfi v N)N, (4. 12) 

J£(A0 =3ff 5 y 5 (y"/ + HvfN 5 - a 4 y 4 fi v (fir + fj v )(l2fj 2 + 35yu /y u v + 12fi 2 v )N 4 

- a 3 y 3 fjl(26fj 4 + I5\fi]fi v + 2Alfi~fi\ + 15 W v + 26fj 4 v )N 3 

- a 2 y 2 fj\{in+fj v ){\2fj\ + %5fj]fj v + I34fijfi 2 v + %5fi,fj 3 v + \2fi 4 v )N 2 

- ayfj 4 v (jj, + fJ V ) 2 (fj 4 + \3fj)fj v + 35fj 2 ,fj 2 v + 13/// fi 3 v + fj 4 v )N - 4fj 2 fj 7 v (fj I + fj v ) 3 . 

(4.13) 

Whereas D(N) is always positive, the sign of "K{N) depends on N. Since 3-£(0) < and 
lirrijv_> +00 Ji(N) = +oo, the function % has at least a positive zero. We define by N t the 
(first) positive zero of J£. Therefore, 3i(N) < for < N < N*. We thus conclude that, for 
Tmax large enough (let us say r max > T^) dx > T^ix, which depends on a, y, fir, fJr and fiy), 
Reci(ri) < for any < N < N*, whereas Reci^) < for any N > 0. This completes 
the proof. □ 



TWO-DIMENSIONAL STABILITY ANALYSIS IN A HIV MODEL 



23 



5. Numerical results 

In order to show the stability of the infected steady state numerically, we can fixiV = 
300, start from the value r e ,.„ o = 0.05625, increase the logistical parameter r monotonically 
until a critical condition is reached such that any further change would result in instability, 
other parameters can be found in Table [1] We present the graphs of numerical solution 
of the system (ll.U -( TT31 l and the trajectory of the solution in the three-dimensional T-V-I 
space. Initial data are To = T u +e(sinxcosy), Io = 0.0, Vo = 0.0185. Some figures assure 
that this solution approaches the limit cycle in the instability subdomain V. In Figure [7] 
corresponding to the subdomain T, the solution approaches the periodic orbit. 





(e) Free virus at (10,10) 



Figure 5. Dynamical solution of System j 1 . 1 b - d 1 .3b on a 20x20 grid. 
Parameter values are N = 300 and r cl ; t = 0.05625 < r = 1 .0 < r\ — 
2.1846. The infected equilibrium is stable. 
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Figure 6. Dynamical solution of System dl.U - dl.3l ) on a 20x20 grid. 
Parameter values are N = 300 and r = 2.0 < r\ — 2.1846. The infected 
equilibrium is stable. 
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Figure 7. Dynamical solution of System dl.U - dl.3l ) on a 20x20 grid. 
Parameter values are N = 300 and r = 200.0 < r 2 = 464.1225. The 
infected equilibrium is unstable. 
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(a) Densities of virus V 




5 10 15 20 25 30 35 
time (days) 

(c) Free virus at (10,10) 





(b) Densities of target cells T 




10 20 30 40 

time (days) 

(d) Target cells at (10,10) 




Virions (V) per mm 3 4 " J Target cells (T) per mm 3 

(e) Free virus at (10,10) 

Figure 8. Dynamical solution of System dl.U - dl.3l ) on a 20x20 grid 
Parameter values are N = 300 and r = 500.0 > r 2 = 464.1225. The 
infected equilibrium is stable. 
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Appendix A. Eigenvalues of the Laplace operator with periodic boundary conditions 

Let A be the realization of the Laplace operator A in l? c , with H 2 C as a domain. The 
following is a well-known result. Nevertheless, for the reader's convenience we provide a 
short proof. 

Theorem A.l. A is a sectorial operator and its spectrum cr(A) is a countable set of 
semisimple eigenvalues. More precisely, 

(T{A) = |~(*? + *!):Jfci,*4eNj. (A.l) 

Proof. Fix A e C, / e l? c and consider the resolvent equation 

Au-Au=f. (A.2) 

Denote by e/, the function defined by 

1 2knt \ 

e h (t) = exp I » t&R, 

for any h 6 Z. Then, the functions (x,y) t-» ef,(x)ek(y) are an orthogonal basis of L^. 
Hence, any function g e Zi, can be expanded into a Fourier series as follows: 



g<X30= Z 1/2 f M ^«3') e -*iW e -fcO'y*rfy)«tiW«jbO') = : u h,ki e kS x )ek 2 (y), 

for almost every (jc, y) e . Multiplying both sides of ( IA.21 > by (x)e,t, (y) and integrating 
over Qf, it thus follows that, if u e fl? c , then the Fourier coefficients of u solves the 
infinitely many equations 



A( 2 u kxM + 4n 2 (k 2 + k 2 2 )u kuk2 = ff kx ,k 2 - 

These equations are uniquely solvable if and only if A + — 75- (^f + k%) and, in this case, we 
have 

fh,h 

Uk,,k 2 



1 + A -fr(k\ + k\) 
A straightforward computation shows that the function 

is in Z)(A) and, actually, solves the resolvent equation, when A + -yr(k 2 + k 2 ,) for any 



k\,k2 e N. We have so proved that <x(A) is given by dA.U . 

It is immediate to check that <x(A) consists of eigenvalues only. Moreover, if Re A > 0, 
we can estimate 

HtfCWII^ = 2 - * 4 Z l^*" 2 = 

Proposition 2.3.1 in [14] implies that A is sectorial in L 2 (Q.(). 

Finally, we show that all the eigenvalues of A are semisimple. For this purpose, let us 
fix one of such eigenvalues Aq and let H = {(fci.fe) e N 2 : Ao = -^-(k 2 + k 2 )}. Then, 

(^-/l )/?(i,A)/(x,y) = ^ fh^h^e^iy) + ^ 4 ^ 2 ^° — -fhte e h( x ) e h(y) 

k u k 2 eH ki,k 2 m ^ + + ^2-* 

=:Pf(x,y) + (A-Ao)Dtf(x,y). 
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Clearly, P is the spectral projection on the eigenspace corresponding to the eigenvalue 

r 2 



/lo. On the other hand, Dx is a bounded operator in Li uniformly with respect to A e 



B(A , 2n l /£ l ). Indeed, if (k x ,k 2 ) $ H and A is as above, then 



4-n 2 

—(k 2 + k 2 2 ) + A 



4-n 2 

—(k 2 +k 2 2 ) + A 



4-n 2 2n 2 
-\A-A Q \> — -\A-A l) \>—. 



Thus, 



i.e., Dx is bounded, uniformly with respect to A e B(Aq, 2jt 2 /{ 2 ). These results imply that 
Ao is a semisimple eigenvalue of A. Note that the eigenspace corresponding to Ao is one- 
dimensional if and only if H is a singleton. In this case, Aq is a simple eigenvalue of A. 
More precisely, the geometric multiplicity of the eigenvalue A = ^-(k 2 + k 2 ,) is given by 
nix = 4 ri/ , =i( r ; where the coefficients r, are given by the following decomposition of 
k 2 + k\ in primes 

m n 

e 1 +%=2«y]p?Y\ q j, 

;=i j=i 

with pi being primes of the form 4t+ 1, and being primes of the form 4f + 3 (see (9 1). □ 

The following classical result on Sturm-Liouville problems is the key tool to prove 
Theoreml2.iriv). 



Corollary A.2. Let d and p be, respectively, a positive constant and a bounded measurable 
function. Further, let B : H 2 —> L 2 be the operator defined by Bu — dAu - fiu for any 
u e H 2 . Then, the spectrum of B consists of eigenvalues only. Moreover, its maximum 
eigenvalue A msx is given by the following formula: 

{ dj m 2 d X + j ^ 2 d X ] 

/W = - mf < J. > . (A.3) 

Finally, the eigenspace corresponding to the eigenvalue A m . ix is one dimensional and con- 
tains functions which do not change sign in Q.(. 
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